Nonlinear resonant behavior of the dispersive readout scheme for a superconducting 

flux qubit 
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the persistent current qubit |21j, |22( , for the charge qubit 
|23| . and for the hybrid qubit where the readout was op- 
erated in the nonlinear regime for its use as a bifurcation 
amplifier [24j. 
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A nonlinear resonant circuit comprising a SQUID magnetometer and a parallel capacitor is 
studied as a readout scheme for a persistent-current (PC) qubit. The flux state of the qubit 
is detected as a change in the Josephson inductance of the SQUID magnetometer, which in 
turn mediates a shift in the resonance frequency of the readout circuit. The nonlinearity and 
resulting hysteresis in the resonant behavior are characterized as a function of the power of both 
the input drive and the associated resonance peak response. Numerical simulations based on a 
phenomenological circuit model are presented which display the features of the observed nonlinearity. 

PACS numbers: 



I. INTRODUCTION 

Superconducting Josephson junction circuits are 
promising candidates for realizing a quantum computer. 
These solid-state qubits can be fabricated using standard 
integrated-circuit techniques, where there is the possibil- 
ity to incorporate the control and readout circuitry on- 
chip, providing a manageable option for scaling up to 
a larger number of qubits. Quantum-coherent phenom- 
ena pj have been studied utilizing the quantum states Q 
of single-qubit circuits and cavities, including superposi- 
tions of distinct macroscopic states IE El , time-dependent 
Rabi oscillations IE Ip , H, I S 0, E cavity quantum- 
electrodynamics G3j Ei 13 and Mach-Zehnder-type in- 
terferometry Coherent oscillations spec- 

troscopic evidence for entanglement 0] , and a prototyp- 
ical gate operation [2(J have also been demonstrated in 
superconducting coupled-qubits. However, to further in- 
crease the coherence times of these qubits for manipula- 
tion of their quantum states, one must find ways to re- 
duce the amount of noise intrinsic to the qubit, as well as 
noise introduced by the readout process itself. In partic- 
ular, several previous readout methods have relied on the 
switching of a Josephson circuit from a zero- voltage to a 
finite-voltage state. This switching generates quasipar- 
ticles, and thus such readout approach is limited by the 
subsequent decoherence. More recently dispersive read- 
out schemes have been developed such that the qubit 
is coupled to a resonator, and the state of the qubit is 
detected as a shift in the resonance frequency of the res- 
onator. As a result, the readout process requires only 
lower input biases and hence minimizes the generation of 
quasiparticles. Furthermore, the resonator also acts as 
a narrow-band filter which shields the qubit from broad- 
band noise. Dispersive readout has been implemented for 
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FIG. 1: Experimental setup: (a) Circuit schematic of the res- 
onant readout circuit. The designed component values were 
Li = 69 nH, L 2 = 0.78 nH, Ci = 1.4 pF, and C 2 = 100 pF. 
The SQUID inductance L j was approximated to be 0.2 nH 
for the circuit design, (b) Optical micrograph of the actual 
device, (c) Electronic setup at different temperature stages 
of the 3 He cryostat. 



This paper focuses on characterizing the nonlinear res- 
onant behavior of the dispersive readout scheme devel- 
oped for a persistent current qubit. The readout ele- 
ment is a SQUID magnetometer, which is operated as a 
nonlinear, flux-sensitive inductor incorporated in an L-C 
resonator. The qubit is coupled to the SQUID inductor, 
and the flux state of the qubit is detected as a shift in the 
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resonance frequency of the resonator by means of magni- 
tude and/or phase measurements. Our approach differs 
from other resonant-type experiments in two main ways. 
First, our qubit and readout circuit were fabricated on 
the same chip from niobium, whereas the implementa- 
tions in |2lj|-|24| were aluminum-based. Second, we were 
able to achieve a high quality factor for the resonator 
by incorporating an RF transformation network on-chip 
using the planarized niobium process. We observed reso- 
nant behavior due to the nonlinear Josephson inductance 
of the SQUID, given the high quality factor of the res- 
onance |25l l2q . The frequency spectra of the readout 
circuit were characterized in both the linear and nonlin- 
ear regimes. Biasing the readout circuit in the nonlinear 
regime potentially provides additional sensitivity for dis- 
tinguishing the qubit states [24| . 

The persistent current (PC) qubit used in this study is 
a superconducting loop interrupted by three Josephson 
junctions, two of which have the same critical current 
while the third junction has a critical current reduced 
by a factor a |27L l2Sj . When the external magnetic flux 
threading the qubit loop is biased near half a flux quan- 
tum, the two lowest energy states correspond to oppo- 
sitely circulating persistent currents in the qubit loop. 
The induced flux of the persistent current (and hence the 
state of the qubit) is detected by a SQUID magnetometer 
which surrounds the qubit. 

In the resonant readout scheme, the SQUID magne- 
tometer is operated using the property that the Joseph- 
son inductance of the SQUID is a nonlinear function of 
both the current bias I sq and the flux bias & ex t- In our 
experiments, the SQUID current bias comprises solely an 
AC component, whereas the flux bias $ e xt = &dc + &ac 
has both a DC component corresponding to the exter- 
nal bias, and an AC component corresponding to the 
induced flux that is mutually coupled to the SQUID. To 
demonstrate the general principles underlying the oper- 
ation, consider the limiting case where the SQUID has 
negligible loop inductance and symmetric junctions each 
with critical current I co . In this limit, the SQUID be- 
haves like an equivalent single junction with an effective 
critical current given by I c {$ext) = 2J co | cos(7r$eat/$ )|, 
and an effective phase given by ip p = sin -1 (I sq j T c (^ext)) ■ 
By defining the Josephson inductance according to V sq = 
d [Ljl sq ] /dt, where V sq is the voltage across the SQUID, 
we obtain the inductance to be 



Lj(I sq ,3> ext ) = 



$ 



27r/ c ($ ext ) sin^p 



(1) 



In the linear regime where the current and flux biases 
are small, one can approximate the SQUID by a linear 
inductor given by 
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(2) 



The inductance for small AC drives can be approximated 
quasi-statically by the inductance at the DC operating 
points for the current and flux biases. To demonstrate 



the separate effects due to the current and the flux, we 
first set ^ext to zero and reduce Eqn. ^ to 



Lj(I sq ,0) = 2L Jo - 
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Lif p 
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where (p p = sm~ 1 (I sq /2I co ). Thus, the inductance in- 
creases with the size of the driving DC current bias. Like- 
wise, when I sq in Eqn. ^ is set to zero, 
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which has a periodic dependence on the DC flux bias with 
periodicity given by <I> . Starting from a bias of Q^c = 0, 
the inductance increases with flux, and starting from a 
bias of $dc = 0.5 $ , the inductance decreases with flux. 

The general trend of an AC bias can be conceptualized 
as averaging the inductance about the DC bias point over 
the range of the AC bias. Hence, near &dc = both the 
AC driving current and the AC flux increase the effec- 
tive inductance as the AC drives increase. In contrast, 
near $,j c = 0.5 $ an increasing AC flux bias tends to de- 
crease the effective inductance and an increasing AC cur- 
rent bias tends to have the opposite effect. Therefore, the 
current and flux act in concert at $>dc = 0; whereas, they 
compete at $dc = 0.5J>n. From our experiments and 
numerical simulations [29j, we have found in our read- 
out circuit that the AC flux dominates at $dc — 0.5 $ G . 
Consequently, this paper will focus on the effects due to 
flux, thereby allowing us to develop a phenomenological 
model that qualitatively matches the experimental ob- 
servations. 

Fig. ^ shows the circuit schematic of the resonant 
readout circuit. The PC qubit is mutually coupled to the 
SQUID inductor Lj. The resonating loop comprises Lj, 
L2, and the parallel combination of C\ and C2. To raise 
the quality factor of the resonance for higher readout sen- 
sitivity, a tapped-inductor transformer formed by L2 and 
Lj is used to step up the effective output resistance at 
the resonance frequency. On the input side, L\ and C\ 
form an L-match network which matches the input re- 
sistance to the transformed output resistance [30|. R s 
and Rl represent the 50 f2 source and load impedances 
from the RF electronics, and no resistors were fabricated 
on-chip. The junctions of the SQUID are each shunted 
by a 5 pF capacitor (not shown). The device was fab- 
ricated using the planarized niobium trilayer process at 
MIT Lincoln Laboratory [3l]]. A device micrograph is 
shown in Fig.^x The Josephson critical current density 
was estimated to be 1.2/iA//mi 2 from the process test 
data. The designed junction dimensions were 1.0 /jm and 
0.9 11m for the qubit, and 1.5 /im for the SQUID. Due to 
process bias, the effective electrical junction dimensions 
are expected to be smaller. We measured the effective 
size of the SQUID junctions to be approximately 1.3 /im, 
a reduction of 0.2 /im from their drawn dimension. The 
effective qubit junction sizes were not measured directly, 
but were estimated to have a reduction of approximately 
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0.35 /im, as determined by measuring similarly drawn 
1.0 /im process-test junctions nearby. The area ratio of 
the SQUID to the qubit loop was designed to be 1.3, 
with mutual coupling estimated to be 30 pH. The induc- 
tors were realized by square spirals with a linewidth and 
spacing of 1 /jm, while the capacitors comprised Nb elec- 
trodes with a dielectric consisting of 50 nm of Nb2C>5 and 
200 nm of Si0 2 . 

Our measurements were taken in a 3 He cryostat at 
300 mK. The measurement setup is shown in Fig.^. The 
DC lines were used to characterize the junction prop- 
erties, while the RF lines were used for the resonant 
readout. The external DC flux bias for the qubit was 
provided by a superconducting coil wrapped around the 
sample housing. The signal from the resonant circuit 
was amplified at room temperature. We measured the 
transmission characteristics of the readout circuit with a 
spectrum analyzer equipped with a tracking generator, 
or with a network analyzer when the phase information 
was needed. We used a resolution bandwidth (RBW) of 
3 kHz, and averaged each spectrum 100 times. 



II. QUBIT READOUT AND EFFECT OF INPUT 
BIAS ON READOUT CIRCUIT 

The resonance frequency of the readout circuit was 
measured to be near 419 MHz, with a quality factor esti- 
mated to be on the order of 1000. Fig. 12k shows the re- 
sults when an external flux bias $dc was applied through 
the sample. At a given <&dc, we measured both the res- 
onance frequency and the peak power of the resonance 
spectrum. The resonance frequency of the readout cir- 
cuit, in the linear regime where all the AC biases are 
small, is related to the effective inductance Lj and ca- 
pacitance C by 



fo(^dc) 



1 



2n^Lj{<S> dc )C 



(5) 



A periodic modulation of the resonance frequency of the 
readout circuit was observed and is interpreted as being 
caused by the periodic modulation of the Josephson in- 
ductance of the SQUID. At every 1.3 times the SQUID 
modulation period, a shift in the resonance frequency, 
corresponding to about 2pH (1%) change in Josephson 
inductance, was observed. These shifts, referred to as 
qubit steps, occur as the flux in the qubit $ g is swept 
past $ g = 0.5 $ , as it is more energetically favorable for 
the qubit to change from one circulating current state to 
the other in order to remain in the ground state. The 
periodicity of the qubit steps (corresponding to a flux 
quantum for the qubit) and the periodicity of the SQUID 
lobes (corresponding to a flux quantum for the SQUID) 
are related by the ratio of their loop areas, which was 
defined by the fabrication parameters. In addition, we 
observed a dip in the resonance-peak power (not shown) , 
which corresponds to a broadening of the resonance near 
the qubit step region [32] • The parabolic-like background 



observed in the frequency modulation curve was due to 
undesired heating from the magnet current in the DC 
(soft-coax) lines. The heating causes an increase in the 
resonance frequency, and is more significant at high mag- 
net current biases. The heating effect was eliminated for 
faster scans and when the sample was later tested in a 
dilution refrigerator using superconducting magnet leads. 
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FIG. 2: (Color online) (a) Modulation of the resonance fre- 
quency with external DC flux bias. Qubit steps are observed 
at 0.18 mA and -0.495 mA. (b) Modulation of the resonance 
frequency for various input power. Amount of modulation is 
reduced in the nonlinear regime (/t.jvL — /o,nl) compared to 
the linear regime (fb.L ~ .fa,L)- The circular markers represent 
the inflection points where d 2 f /d& 2 = 0. 



Fig. |2t> shows that as the level of input bias increases, 
the amount by which the frequency is modulated over a 
flux quantum decreases. This will be shown in the next 
section to be a direct consequence of the shape of the 
resonance spectrum as it becomes increasingly nonlinear 
with higher input power. 



III. NONLINEAR RESONANT BEHAVIOR OF 
READOUT CIRCUIT 

The resonant readout circuit can experimentally dis- 
tinguish the difference in the flux produced by the circu- 
lating current states of the qubit. Given that the efficacy 
of the readout scheme depends on the nonlinear response 
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of the readout circuit, we now characterize the resonant 
behavior of the readout circuit as a function of the DC 
flux bias for higher AC drives. 

Fig-Blshows the evolution of the magnitude and phase 
spectra with increasing input power for external flux bi- 
ases of $d c = 0, 0.3 $ and 0.5 <Jv In the case when 
&dc = 0, the magnitude and phase spectra evolve from a 
symmetric shape to being asymmetric with a lower res- 
onance frequency as the power of the drive is increased. 
The lower resonance frequency indicates that the effec- 
tive Josephson inductance over an oscillating period is 
higher. For higher levels of the input power, the mag- 
nitude spectrum exhibits a discontinuity near the reso- 
nance frequency, where the system jumps from the lower 
branch to the higher branch. The phase spectrum also 
exhibits a discontinuity similar to the magnitude spec- 
trum. For $dc = 0.5 $ , the asymmetry is opposite to 
that of &dc = 0; the resonance frequency increases with 
higher power, indicating that the overall effective induc- 
tance is decreasing with increasing power of the drive. 
An intermediate behavior is captured at $dc = 0.3 $ . As 
the input power increases, the nonlinear magnitude spec- 
trum first bends towards the lower frequency side, then 
gradually evolves into a characteristic shape with two dis- 
continuities near the resonance frequency, once when the 
magnitude is increasing and once when the magnitude 
is decreasing. Similarly, the phase spectrum also shows 
two discontinuities at the same frequency locations, with 
a partial phase drop at each discontinuity. 




Frequency [MHz] Frequency [MHz] Frequency [MHz] 

FIG. 3: Evolution of the magnitude and phase spectra of the 
readout circuit from the linear to the nonlinear regime with 
increasing input power. Data are shown for flux biases at 
&dc = 0, 0.3 <3? and 0.5 $ . The nonlinear spectrum evolves 
from having a lower resonance frequency at $ = to having 
a higher resonance frequency at $ = 0.5 <& D . A self-resonance 
(due to parasitic couplings) was observed near the resonance 
frequency of the spectrum. The phase spectrum at -54dBm 
for $ = 0.3 <I? was arbitrarily shifted for display purpose. 

The shapes of these curves are similar to the response 



of driven, weakly nonlinear systems which exhibit an in- 
stability region indicating multiple solutions and hystere- 
sis (33[[34|. Two such curves are shown in Fig.0Ji and b. 
In particular, we model our system as a nonlinear circuit 
which results from a current-driven LRC resonant circuit 
with a nonlinear inductor L. In this case the flux in the 
inductor <f> satisfies 

d 2 $ 1 ri$ 
Ismw s t = C-^ + -— +h($,d$/dt), (6) 

where the function h(Q,d$/dt) models the nonlinearity 
of the inductor. For example, when h = &/L a then 
the system is a simple LRC resonant circuit with a lin- 
ear inductor L . When h ~ sin<I> the nonlinear equa- 
tion is analogous to a driven pendulum system whose 
response is similar to Fig. 0J) |33l l37| . Another exam- 
ple is the Duffing Equation where h ~ $ — c$ 3 , whose 
response is like Fig. J2Ji for negative c and like Fig. 0Jd 
for positive c l34|. In section [V] we will use a func- 
tional form for the effective inductance which incorpo- 
rates both the needed dependence on applied DC flux and 
the resonance-frequency dependence observed for small 
drives. In fact, given that Fig. shows that the reso- 
nance frequency is periodic in the applied DC flux, then 
the effective inductance that needs to be captured in the 
form of h must also follow this same periodicity. 
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(a)o dc =o.5<s o (b)4> dc = o 




Frequency 



1 1 : Frequency 

f a,L f b,L 

FIG. 4: Illustration of the resonance spectra in the nonlin- 
ear regime (a and b) and the linear regime (c and d) for 
$ = 0.5 $0 and <3> = respectively. The shaded region of 
the nonlinear spectrum marks the region over which multiple 
solutions occur. The solid line traces the actual spectrum ob- 
served experimentally with a forward frequency sweep, and 
the circular marker corresponds to the peak frequency that 
was being measured. The bending of the nonlinear spectra 
to opposite sides accounts for the reduced separation of reso- 
nance frequencies (fb,NL — .fa,NL) compared to the linear case 
(fb,L — /o,i), as was observed in Fig. |5J>. 

We now use the illustration in Fig.0]to explain some of 
the general features of the data in Fig. [21 and to motivate 
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the subsequent analysis. In Fig.QJ the shape of the reso- 
nance spectra are shown for both the linear and nonlinear 
regimes for $dc = and 0.5 <I> , with the resonance fre- 
quency at $dc — 0.5 $ lower than at &dc = given the 
flux-dependence of the effective inductance. The shaded 
region of the nonlinear spectrum corresponds to the re- 
gion / S [fup, fbend] over which multiple solutions occur 
(two of which are stable and one of which is unstable) 
|35| . The solid line traces the actual spectrum observed 
experimentally with a forward frequency sweep, and the 
circular marker corresponds to the peak frequency that 
was being measured. 

We have seen from Fig. [5p that as the level of input 
power increases, the amount by which the frequency is 
modulated over a flux quantum decreases. This is a di- 
rect consequence of the shape of the resonance spectrum 
as the system response becomes increasingly nonlinear. 
As illustrated in Figs. and d, the resonance spectra 
at <i> = and 0.5 $ have resonance frequencies that are 
maximally separated (fb.L ~ fa,i) when the input bias is 
low, and therefore when the resonance spectra are nearly 
those of a linear response. As the input current bias 
increases, the resonance spectrum evolves from the sym- 
metric Lorentzian shape to an asymmetric shape. This 
is shown in Figs. 0Ji and b. The fact that the nonlin- 
ear spectra bend to opposite sides at $ = and 0.5 $ 
accounts for a reduced amount of modulation in reso- 
nance frequency (fb,NL — fa,Ni) compared to the linear 
case. It was also observed in Fig. \Bp that the frequency 
modulation curves for different input power meet period- 
ically at the inflection points, where the second deriva- 
tive d 2 f /d§ 2 equals zero. In fact, the asymmetry of the 
spectrum changes sign near the inflection points. 

To further quantify the amount of bending in the non- 
linear spectrum, we introduce a parameter Sf which is 
a normalized shift of the resonance frequency /„ of the 
nonlinear spectrum relative to the linear spectrum f a : 



Sf = 



In fo 



(7) 



Experimentally, f a was determined as the resonance fre- 
quency of the spectrum measured at the lowest power 
(-74dBm). /„ was defined as the peak frequency and, 
in the limit of high input power, the frequency at which 
the discontinuity occurs. The sign of Sf serves as an in- 
dication of the polarity of the bending. A positive Sf 
corresponds to the nonlinear spectrum bending to the 
higher frequency side, and a negative Sf corresponds to 
the spectrum bending to the lower frequency side. 

In Fig.[S]the normalized frequency Sf of the resonance 
spectrum is plotted for increasing input power from - 
74dBm to -54dBm. The measurements of the spectra 
were made with a forward frequency sweep. The differ- 
ent markers correspond to various flux biases between 
$ dc = to 0.5 $ . At $ dc = 0.5 $ (top plot), Sf 
is increasingly positive; whereas, at $dc = (bottom 
plot), Sf becomes increasingly negative. Furthermore, 
the amount of bending \Sf\ at $>dc = is smaller than 



at 0.5 $ for a given input bias, which is related to the 
fact that a forward frequency sweep captures the full fre- 
quency extent of the bistable region for &dc = 0.5<I> o 
but not for <&dc = 0. (The reverse is true if the fre- 
quency is swept backwards, as discussed in the next sec- 
tion.) Finally at intermediate flux biases between 0.3 <f> 
to 0.36 $0, Sf shows an undulating behavior, correspond- 
ing to the asymmetric spectrum constantly varying its 
polarity. There is an initial linear dependence of Sf on 
input power, which will be discussed in the next section. 
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FIG. 5: (Color online) Bending Sf as a function of input 
power from -74dBm to -54dBm in ldB steps. Measurements 
were made with a forward frequency sweep. The data are 
shown for various DC flux biases between to 0.5 $ . The sign 
of 5f indicates the polarity of the asymmetric spectrum. The 
red line corresponds to a linear fit in the low-power regime. 



IV. HYSTERESIS OF THE RESONANCE 
SPECTRUM 

For larger input drives, the resonance spectrum ex- 
hibits a discontinuity which corresponds to one of the two 
boundaries of the bistable region. Within the bistable re- 
gion, the system settles into one of the solutions depend- 
ing upon the initial conditions. For our case, the initial 
condition is set by the solution at the previous driving 
frequency, which in turn is determined by the direction 
of the frequency sweep. The resonant behavior of the 
readout circuit presented so far were obtained with a for- 
ward frequency sweep. Here, we present the hysteretic 
behavior of the resonance spectrum measured with both 
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forward and backward frequency sweeps so that the full 
boundary of the bistable region can be mapped. 
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FIG. 6: (Color online) Top plots: Hysteretic resonance spec- 
trum for flux biases at <&dc = 0, 0.3 $ and 0.5 $ D . The extent 
of the bistable region is given by |/t — ft\- Middle plots: 8 ft 
and 8 ft as a function of input power from -66 dBm to -54 dBm 
in ldB steps. Onset of hysteretic regime occurs at Ph=- 
61 dBm for $ dc = 0.5 $ D and P H =-59dBm for $ dc = 0. Bot- 
tom plots: 8 ft and 8ft re-plotted as a function of resonance- 
peak power. The dotted line is a linear fit for low power. 



The top three plots in Fig. show the typical hys- 
teretic spectrum for flux biases at $dc = 0, 0.3 $ and 
0.5 $ . The data are shown for an input power level of - 
54 dBm, which corresponds to a highly nonlinear regime. 
The direction of the frequency sweep is indicated by the 
arrows. For the case of $dc — and 0.5 $ , we define 
the extent of the bistable region as \f t — fb\, where fb 
is the frequency at which the resonance spectrum jumps 
from the lower to the higher stable branch, and ft corre- 
sponds to the frequency at which the spectrum falls from 
the higher to the lower stable branch. At $dc — 0.3 $ , 
the forward and backward traces overlapped, indicating 
that the bistable region associated with the two discon- 
tinuous edges were too small to be detected given the 
frequency resolution. 

The onset of the hysteretic regime is illustrated in the 
middle three plots of Fig. [U where the extent of the 
bistable region was characterized as a function of input 
power from -66 dBm to -54 dBm in 1 dB steps. First, we 
normalized ft and fb with respect to the resonance fre- 
quency f in the linear regime according to a definition 



similar to Eqn. [7] 

ft — fo 



5ft 



fo 



and Sf b = ^-A (8) 

fo 



5ft and Sfb were then plotted as a function of input 
power. Hysteretic behavior was observed when the in- 
put bias was above a threshold Pjj, which was measured 
to be -61 dBm for $ = 0.5 $ Q , and at a higher power of 
-59 dBm for $ = 0. 

In the bottom row of plots of Fig. El we have plotted 5ft 
and 5fb as a function of the resonance-peak power. We 
see that 5ft is initially a linear function of the resonance- 
peak power for all three fluxes shown and that it is fully 
linear for zero DC flux. This initial linear dependence on 
resonance peak power is found for many functional forms 
of the nonlinear term h(&,d&/dt) [33,[34|]. However, the 
dependence on input power varies according to the par- 
ticular functional form of h(&,d&/dt), see reference [34| 
for some sample cases. For example, both the Duffing 
equation and the pendulum model give a linear depen- 
dence on both input and resonance-peak powers in the 
weakly nonlinear regime 0HHH1]. For h = &(d$/dt) 2 , 
5ft depends linearly on the resonance-peak amplitude, 
but as the cube root of the input power [34|. 



V. SIMULATIONS OF NONLINEAR 
RESONANT BEHAVIOR BASED ON A 
PHENOMENOLOGICAL MODEL 

In this section, we present simulations to illustrate the 
nonlinear resonant behavior of the readout circuit. One 
approach to simulate the AC-driven behavior of a cir- 
cuit comprising a SQUID is to numerically solve the set 
of coupled differential equations governing the SQUID 
consistently with the rest of the circuit. However, the 
dynamical modeling of the resulting circuit is complex; 
for example, for a circuit with a SQUID shunted by a 
resonating capacitor has 6 dynamical variables when the 
mutual inductive coupling between the SQUID and the 
resonating loop is included [2^|. Therefore, we use the 
phenomenological LRC circuit model of Eqn. with the 
linear inductance replaced by a flux-dependent nonlin- 
ear inductor. This approach allows a reduction of the 
mathematical complexity of the problem to one dynam- 
ical variable. 

Specifically, the readout SQUID is modeled by a flux- 
dependent nonlinear inductor Lj(& ext ) given by 



Lj{<5>ext) 



Lo 



y/(l + P 2 ) + (1 - j3 2 ) cos(2^$ e:ct /$ ) 



(9) 



The functional form for the nonlinear inductor Lj(<fr ex t) 
in Eqn. [5| captures the sinusoidal-like shape of the fre- 
quency response of the actual readout circuit as previ- 
ously shown in Fig. Eh- This is also illustrated in Fig. [3d, 
where Lj($> ext ) is plotted for L = \/2Lj Q (Eqn. |2J, and 
for different values of (3. It can be seen that (3 has an 
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is then governed by Eqn. H3 and 



-- p 


= 


p 


= 0.3 


— p 


= 0.5 



1 ^extt'n^o] 



FIG. 7: (a) Circuit schematic of the phenomenological res- 
onant model, (b) A plot of Lj(& ext ) given by Eqn. H for 
different values of j3. The circular markers represent the in- 
flection points where d 2 Lj/d& 2 is zero. The illustration shows 
that depending on the DC flux bias, the AC modulation of 
Lj due to $ nc can result in a lower (&dc = 0.5 4> ) or higher 
(&dc = 0) effective inductance . 



effect on (a) the amount by which Lj is modulated over 
half a flux quantum, and (b) the position of the inflec- 
tion points at which the second derivative d 2 Lj/d& 2 is 
zero. One can therefore use j3 in Eqn. as a fitting pa- 
rameter such that the location of the inflection points of 
the frequency response (d 2 / /d$ 2 = 0) match the data 
in Fig. ^jp. The functional form of the inductance was 
motivated by that of an asymmetric SQUID j3?J since it 
has the needed periodicity; however, it is only the form 
that is used, and we are not assuming that asymmetry 
plays any role in the physical circuit. 

We further assume that the inductor is mutually cou- 
pled to a total external flux bias as Q e xt = &dc + &ac, 
where <&dc is the DC flux bias that was applied experi- 
mentally to the SQUID, and $ ac is any oscillating flux 
that was mutually coupled to the SQUID. This AC-part 
of the flux is modeled as being proportional to the self- 
induced flux as $ oc = u<& to ensure that the amount of 
coupled flux increases as $ gets large near the resonance 
frequency of the circuit. The size of <i> ac used in the sim- 
ulations is around 0.1 <&„ near the resonance frequency. 

The dynamics of the phenomenological resonant circuit 



(10) 



For small drives (a m 0), this term can be approx- 
imated by a linear inductance which depends on <&dc 
in a sinusoidal-like fashion. For larger values of the 
drive, Eqn. 1101 can be expanded to be in the form of 
a$ + 6 $ 2 + c$ 3 when (1 - /3 2 )/(l + /3 2 ) << 1. More- 
over, b ~ sin (27r$ dc /$ ) and c ~ cos (27r$ dc /$ ). For 
$>dc = and $dc = 0.5 $ , then 6 = 0, and the resulting 
equation is of the form of a Duffing Equation with both a 
linear and cubic term in so that with increasing drive 
there will be a bending of the resonance frequency with 
its associated hysteresis [34[ ■ Moreover, the sign of the 
cubic term is opposite for $dc = and $dc — 0.5 $ Q , and 
hence bending will be in opposite directions, as needed 
to qualitatively explain the data. At other values of $dc 
there will be a competition between the quadratic and 
cubic terms. 
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1.05 0.95 
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FIG. 8: (Color online) Simulated magnitude and phase spec- 
tra of Vl for increasing drive amplitude /. The results qual- 
itatively reproduce the experimentally observed behavior in 
Fig. |3] Vl is plotted in reduced units of <& / \/ Lj a C and 7 
in units of 21 co . The frequency axes are normalized with re- 
spect to the resonance frequency of the linear spectrum: f a — 
1.88 GHz ($ dc = 0), 1.59 GHz (0.3 $ D ) and 1.33 GHz (0.5 $„). 
The phase spectra at the highest drives I = 0.1 and 0.2 for 
$dc = 0.3 $ are arbitrarily shifted for display purpose. 



By numerically solving Eqn. for $ at different driving 
frequencies oj s , the magnitude and phase spectra of the 
voltage across the inductor Vl = d<fr/dt were obtained. 
In Fig. [8] the spectra are shown for increasing drive am- 
plitude for <J> dc = 0, 0. 3 <& a nd 0.5 $ D . Vl is plotted in 
reduced units of & /^L Jo C, and the drive amplitude I 
is in units of 2I co . The driving frequency was swept such 
that the lower stable branch within the bistable region 
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is shown for all flux biases. The simulation qualitatively 
resembles the experimental data presented in Fig. As 
expected, the nonlinear spectrum has a lower resonance 
frequency at (&dc — 0, and a higher resonance frequency 
at = 0.5 <£> . Also, discontinuities are observed at the 
boundary of the bistable region for the higher biases. At 
&dc = 0.3 $ , the magnitude spectra at the two highest 
input biases exhibit two discrete jumps, once at a lower 
frequency when the magnitude is increasing, and another 
at a higher frequency when the magnitude is decreasing. 
As for the phase spectra, a partial phase drop occurs at 
the low-frequency discontinuity while most of the phase 
drop occurs at the high-frequency discontinuity. 

We have also performed simulations which qualita- 
tively reproduce the hysteresis data in Fig. This was 
done by stepping the driving frequency in both the low- 
to-high and high-to-low frequency directions, and by en- 
suring the initial conditions used for the next frequency 
point were the solutions obtained for the previous fre- 
quency point. We have assumed that the square of the 
drive amplitude I 2 for the simulations is proportional to 
the input power for the experiment. 

The top row of plots in Fig. El shows the typical 
simulated hysteretic behavior at various flux biases for 
I/2I co = 0.01. The extent of the bistable region given 
by \Sfb-Sft\ is the largest at 3>dc = 0.5 $ for this drive 
amplitude. The middle and bottom rows of plots show 
the dependence of Sfb and 8ft as a function of input 
power and resonance-peak power respectively. The trend 
at $d c = 0.3 $ shows qualitative resemblance to the ex- 
perimental data. The magnitude of Sf t initially increases 
linearly with the power of the drive for $dc = and 0.5 $ 
similar to the data. As explained previously, this linear 
dependence is expected at low drives due to the nature 
of the nonlinearity 0, ■ 

The phenomenological model presented here is meant 
to show the qualitative trends in the data. To be more 
quantitative, we have analyzed more complex circuits 
|29| . For example, we have considered the current- 
driven circuit across an asymmetric SQUID with self- 
inductance, and mutual inductive coupling between the 
SQUID loop and the resonating loop. The induced flux 
in the resonating loop is enhanced by the quality factor 
Q near the resonance frequency. The simulation gives 
results which qualitatively reproduce the data with rea- 
sonable numbers; however, the quantitative fitting of the 
data was not possible due to the uncertainty in the actual 
on-chip values for the capacitances and the mutual/self 
inductances. 



VI. DISCUSSION 

In this paper, we characterized the resonant behavior 
of the readout circuit to be utilized in a resonant scheme 
for measuring a PC qubit. We have identified the lin- 
ear, nonlinear (bistable), and hysteretic (highly nonlin- 
ear) bias regimes of operation. Given the high quality 
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FIG. 9: (Color online) Top row: simulated hysteretic behavior 
at various flux bias for drive amplitude I = 0.01. Middle row: 
5 ft and 8 ft as a function of the square of the drive amplitude, 
which is proportional to the input power. The simulations 
were performed for I between 0.001 to 0.05. Bottom row: Sfb 
and Sft as a function of the square of the voltage response 
Vl, which is proportional to the resonance-peak power. Vl is 
in reduced units of 5? /V Lj C and / is in units of 2I co . 



factor of the resonance, we observed a manifestation of 
the nonlinearity due to the Joscphson inductance of the 
readout SQUID. The resonance spectrum of the readout 
circuit became asymmetric in the nonlinear regime, and 
the polarity of the asymmetry changed sign as a function 
of DC magnetic flux bias to the SQUID. The numeri- 
cal simulations using a nonlinear inductor qualitatively 
reproduce the trends in the experimental data. 

To perform time-resolved measurements of the qubit 
on a microsecond timescale, the resonant readout is to 
be operated at a bias frequency f s near the resonance 
frequency f . The change in the resonance frequency 
due to the qubit signal is thus detected as a difference in 
the magnitude or phase of the output voltage at f s . 

Finally, it should be mentioned that while operating 
the resonant readout in the linear regime keeps the in- 
put bias low and reduces the level of decoherence on the 
qubit, the readout operated in the nonlinear regime has 
the advantage of being used as a bifurcation amplifier 
|24| . In particular, the bias frequency f s can be chosen 
within the bistable region of the nonlinear spectrum such 
that the system has two stable solutions corresponding 
to different voltages. The probability of occupancy in 
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the higher/lower stable solution is sensitive to changes 
in the resonance frequency f a (qubit-mediated) relative 
to f s , and thus provides additional sensitivity for qubit 
readout over the linear approach. 
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